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Active (catalysed) and passive (intrinsic) nucleosome repositioning is known to be a crucial event 
during the transcriptional activation of certain eucaryotic genes. Here we consider theoretically 
the intrinsic mechanism and study in detail the energetics and dynamics of DNA-loop-mediated 
nucleosome repositioning, as previously proposed by Schiessel et al. (H. Schiessel, J. Widom, R. F. 
Bruinsma, and W. M. Gelbart. 2001. Phys. Rev. Lett. 86:4414-4417). The surprising outcome 
of the present study is the inherent nonlocality of nucleosome motion within this model - being 
a direct physical consequence of the loop mechanism. On long enough DNA templates the longer 
jumps dominate over the previously predicted local motion, a fact that contrasts simple diffusive 
mechanisms considered before. The possible experimental outcome resulting from the considered 
mechanism is predicted, discussed and compared to existing experimental findings. 

PACS numbers: 



I. INTRODUCTION 

The nucleosome, the most abundant DNA-protein 
complex in nature, is the basic unit of eucaryotic chro- 
matin organization. It is roughly a cylinder of 6nm height 
and lOnm diameter, consisting of a protein octamer core 
and 147 basepairs (bp) of DNA tightly wrapped around 
it in 1 and 3/4 left-handed super helical turns. The genes 
of all higher organisms, ranging from simple ones like 
yeast to most elaborate like humans are all organized 
in long arrays of nucleosomes with short DNA segments 
(linkers) of 50-100 bp interpolating between them, com- 
parable to a beads-on-a-string chain ^ . The higher 
order organization of these units, being most probably a 
solenoid- or zig-zag, crossed-linker- like fiber with 30 nm 
diameter is still under great dispute though it received in- 
creasing theoretical and experimental support in recent 
years. Above that scale of organization, the higher or- 
der structures which link the 30 nm to the final "big X" 
like structure, the packed chromosome, are still unknown. 
Though there are several biologically motivated specula- 
tions about the "big X" its definite structure remained a 
long lasting puzzle for the last 20 years, defying all bio- 
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FIG. 1: The basic problem setting: how does the histone- 
octamer move along the DNA template? Below: the DNA 
loop mechanism as proposed in in Ref. [Il0| 



physical, biochemical and molecular genetics efforts to 
resolve it because of its intrinsic softness and fuzziness. 

An additional obstacle for understanding the chro- 
matin structure is the fact that it is highly dynamic on 
all organization scales. Starting at the macroscopic chro- 
mosome level we see that its structure can strongly vary 
throughout the cell cycle on timescales of hours or days. 
Below that on timescales of seconds and minutes, the 
structure of the 30 nm fiber itself is subjected to great 
variations due to transcription, replication, biochemical 
modification and other dynamic processes. Finally, at the 
lowest organization level, the nucleosome itself has been 
shown to be a dynamical structure being moved along the 
DNA by chromatin remodeling complexes on expense of 
ATP H,Jql . Interestingly, it was experimentally observed 
lH' 0i if ^that nucleosomes can move even autonomously 
on short DNA segments. This intrinsic repositioning be- 
havior was shown to be strongly temperature dependent. 
At room temperature it occurs roughly on timescales of 
~1 hour indicating the existence of significant energetic 
barriers. Besides the fact that the repositioning does 
indeed occur and is of intramolecular nature (the nucle- 
osome stays on the same DNA segment) the underlying 
scenario could not be figured out. It was speculated by 
Pennings et. al. j|, ^ that the mechanism was some 
type of nucleosome-sliding or screwing motion. An alter- 
native explanation which appears to be more consistent 
with the discrete jumps and large barriers observed by 
Pennings et. al. has been recently proposed in Ref. pO[ . 
In this model the basic step in the repositioning process is 
a partial unwrapping of DNA from the very ends of the 
nucleosome 
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|12| followed by a backfolding of DNA 
with a small 10 bp mismatch (cf. Fig. 1). The result of 
this process is the formation of a small DNA bulge or loop 
on the octamer surface. Once trapped on the nucleosome 
surface this small defect carrying some discrete quantum 
of DNA extra length (a multiple of 10 bp, the DNA re- 
peat length) can propagate by diffusion in both direc- 
tions. If the loop happens to surround the nucleosome 
and comes out at the opposite side (in respect to where it 
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was created) the nucleosome is eventually repositioned by 
a distance given by the "pulled in" extra length. The en- 
ergetic barrier and rates of repositioning were computed 
and were shown to be consistent with the Pennings et 
al. experiment [|[ ^. Moreover, the 10 bp discrete step 
repositioning observed in the experiment (discrete bands, 
no 1 bp spaced intermediates) came out as a natural con- 
sequence of the loop length quantization. The latter is 
enforced by the strongly preferred DNA minor groove - 
octamer interaction and the discrete binding sites at the 
nucleosome surface as deduced from the crystallographic 
structures |14| . 

In Ref. |L0[ small loops with short excess length of typ- 
ically — 2x10 bp were considered and it was shown 
that the looping energies involved increase rapidly with 
the excess length implying that only the shortest (10 bp) 
loop contributes significantly to the repositioning mech- 
anism. Consequently the model predicts a classical dis- 
crete random walk with a jump-size of 10 bp instead 
of a 1 bp motion that would be implied by sliding/cork- 
screwing mechanism. Apart from the discrepancy in the 
elementary step size, both models predict very similar 
behavior: a local one-dimensional diffusive motion along 
the DNA chain. 

In this paper we will carefully reanalyze the idea of 
loop-mediated repositioning by applying the classical tool 
of the Kirchhoff kinetic analogy which provides us with 
analytic solutions of the loop problem and enables us to 
look at loops of virtually any given excess length. The 
main outcome of our study will be a different picture of 
repositioning which physically results from the looping 
mechanism: on short up to moderately long segments of 
up to 2-3 X Ip [Ip: DNA persistence length) the repo- 
sitioning is a jumpy process with largest possible loops 
being the most dominant ones in contrast to short 10 bp 
steps as conjectured before. For longer and very long 
(infinite) DNA segments there is an optimal jump size of 
order ^ 0{lp) and the behavior is superdiffusive in con- 
trast to the previously predicted diffusive mechanism. As 
we will see below, these predictions allow us to clearly 
distinguish between different repositioning mechanisms 
in experiments expected to be performed in near future 

II. ENERGETICS OF LOOPS 

Let us now consider the energetics of an intranucleoso- 
mal DNA loop. We will describe it within the framework 
of the Euler-Kirchhoff theory for the static equilibrium of 
rods (Fig. 2). For simplicity and because of the approxi- 
mate planarity of the problem we can in first approxima- 
tion assume the nucleosome and the loop-forming DNA 
to be in one plane and the DNA to be free of any twisting 
deformation. It this case the free energy of our system 
is simply divided into two components, the planar elastic 
DNA-bending and a histone-octamer DNA interaction: 

Utot — Ubend + Uads (1) 




FIG. 2: The Kirchhoff kinetic analogy between the spinning 
top and the bent/twisted rod depicted for a special case: the 
plane pendulum - planar rod equivalence. The inset shows 
how an intranucleosomal loop can be constructed by inscrib- 
ing the octamer (gray disk) into the bent rod. The nucleosome 
opening angle 2a accounts for the adsorption energy cost (see 
text for details). 

The bending energy (within the linear elasticity approx- 
imation) can be written in terms of the local DNA cur- 
vature K 

Ube7id = 77/ (s) ds (2) 

^ J-L/2 

with A w 50 mm ■ ksT being the bending rigidity of 
DNA at room temperature and physiological salt con- 
centrations . The rod is assumed to be parametrized 
by its contour length parameter s ranging from —L/2 to 
L/2 with L being the total length of the loop. The latter 
can be expressed in terms of two independent quantities: 
the excess length AL and the nucleosome opening angle 
a (Fig.2) 

L{a,AL) ^2aR + AL (3) 

where i? « 4 nm is the effective nucleosome radius, or 
more precisely the distance from the center of the nucle- 
osome to the central DNA axis. Because the DNA can 
enter the nucleosome only in quantized orientations (with 
its minor groove phosphates) and bind only to discrete 
positions on the protein surface [Q, the excess length 
AL — n X hp)NA is to a good approximation an integer 
multiple of the DNA repeat length hptNA — iAnra. 

The second part in the total energy Eq. Uads comes 
from the (predominantly electrostatic) interaction be- 
tween the positively charged protein surface and the neg- 
atively charged DNA. It can be roughly measured from 
experiments probing the competitive protein binding to 
nucleosomal DNA |jri|, |l2j . Neglecting the discreteness of 
charges (binding sites) on the histone octamer surface it 
can in first approximation be assumed to be proportional 
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to the opening angle a and the adsorption energy density 



ds 



(4) 



with Eads ~ 0-5 — 1-0 kBT/nm as roughly extracted from 
^. Here and in the following we assume an interme- 
diate value of Eads = 0.7 kBT/nm. 



A. Ground states of trapped loops 

In order to compute the ground state for a trapped 
intranucleosomal loop we have to consider shapes that 
minimize the total energy^ under two constraints: 

1. The excess length AL is prescribed. Therefore we 
have the relation Eq. |^ between the opening angle 
and the total loop length L 



AL = L — 2aR = const. 



(5) 



2. At the two ends s± = zLL/2 the rod has to be tan- 
gential on an inscribed circle of given radius (rep- 
resenting the nucleosome) ^: 



R = 



— const. 



(6) 



Here x (s) and y (s) are the Cartesian coordinates of 
the rod axis as a function of the arc-length parameter s 
(cf. Fig. 2). The absolute value in the second constraint 
needs to be introduced formally for dealing with crossed 
rod solutions (which we consider later on) and can be 
omitted for simple uncrossed loops. 

For an analytical description it is convenient to use the 
angle 9 = 9 {s) between the DNA tangent and the y axis 
(cf. Fig. 2) as a variable describing the DNA centerline. 
In this case the integrated sine (cosine) of 9 over the arc- 
length parameter s gives the x (y) Cartesian coordinate of 
any point along the rod, and the squared derivative (9') 
gives the rod curvature k. Furthermore the nucleosome 
opening angle a is simply related to 9 at the boundary 



9 {L/2) for simple loops 
TT — 9 {L/2) for crossed loops 



The two constraints Eq. ^ and Eq. ^ can be rewritten in 
terms of 9 and then be introduced into the minimization 



by two Lagrange multipliers ^J-l/2■ We then arrive at the 
following lengthy functional 



= A 



L/2 



{9'fds +2aReads 



-m [L - {AL + 2aR)] 

f.L/2 



+A*2 



COS 9ds — R sin a 



(7) 



Here the first line is the bending + adsorption energy 
contribution, the second and third line are the imposed 
length and tangency constraint. Eq.|^ can be rearranged 
in a more familiar form 



L/2 



/i2 COS 9j ds + b.t. 



(8) 



Here b.t. denotes the boundary terms (depending on 
9 {L/2) only) that obviously do not contribute to the first 
variation inside the relevant s interval. The integral in 
Eq.^ is evidently analogous to the action integral of the 
plane pendulum with A {9')^ corresponding to the kinetic 
and —^2 cos 9 to the potential energy of the pendulum. 
The latter analogy is a rather unspectacular observation 
knowing the celebrated KirchhofF's kinetic mapping be- 
tween deformed rods and the spinning top, which con- 
tains our present problem as a simple special case. The 
Kirchhoff's analogy states that the equilibrium confor- 
mations of (weakly) deformed thin rods can be mapped 
to the time-dynamics of a heavy symmetric spinning top 
subjected to a gravitational force. It has been repeatedly 
applied (with or without direct reference to Kirchhoff ) to 
DNA related problems during the last 20 years (e.g. see 
|16|, 0, |l|, 1^, H, 
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23|, |2J). For a nice visual re- 
view on the spinning top-elastic rod analogy the reader 
is referred to Ref. |^ where the general solutions to- 
gether with a "kinetic dictionary" (time t < — > length 
parameter s, gravitational force< — s-rod tension /127 axis 
of revolution< — >tangent vector etc.) are also provided. 

The nice thing about Kirchhoff's analogy apart from 
its esthetical content is that it provides us with explicit 
expressions for DNA shapes subjected to twist, bend- 
ing and various geometric / topological constraints. In 
our simple planar and twistless case, the "spinning top" 
simply reduces to the simple plane pendulum. The corre- 
sponding planar and twistless rods, also called the Euler 
elastica, are most generally given by 



cos^ (s) = 1 — 2msn^ ^— | nij 



(9) 



which can be integrated to obtain the general planar rod 
shape in Cartesian coordinates: 



^ In Eq. ^ we assume that the interaction is only short ranged 
(contact interaction) which is justified by the very short Debye 
screening length of Inm under physiological salt conditions. 

^ Because of the symmetry we have to impose the conditions only 
on one side. 
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(10) 
(11) 
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with sn, dn, cn(. | m) being the Jacobi elliptic functions 
with the parameter m and 



E (u I m) 



dn^ (v I m) dv 



(12) 



denoting the incomplete elliptic integral of the second 
kind in its "practical" form''. The two parameters m > 
and A > in Eqs. |l^, ^ characterize the shape and 
the scale of the solution, respectively. These solutions 
are up to trivial plane rotations, translations, reflections 
and shifting of the contour parameter s ^ s + sq the 
most general planar Euler elastica corresponding to the 
plane pendulum. For different parameters m one obtains 
different rod shapes corresponding to different solutions 
of the spinning top (plane pendulum) motion p5| . The 
case m = describes a pendulum at rest correspond- 
ing to a straight rod, for < m < 1 one has strictly 
oscillating pendulums corresponding to point symmetric 
rod shapes which for m < 0. 92 are free of self inter- 
sections like the one depicted in Fig. 2. For m higher 
than 0. 92 the rods show varying complexity with a mul- 
titude of self-intersections and for m = 1 one has the 
so-called homoclinic pendulum orbit corresponding to a 
rod solution having only one self intersection and becom- 
ing asymptotically straight for s ±oo (for details see 
Ref. [^). For even higher values * of m, i.e., for m > 1 
we have revolving pendulum orbits corresponding to rods 
with self-intersections lacking point symmetry. Finally, 
the limiting case m ~* oo corresponds to the circular rod 
shape. 

In order to describe a trapped loop we need to use Eqs. 
and |l^ imposing the constraints Eq. || and Eq. ||. 
It turns out to be more convenient to replace the pa- 
rameter set (A, m, L) with the new (but equivalent) set 
{X,m,a :— j^) where we introduced the new dimension- 
less parameter a which we call the "contact parameter"^. 
From Eq. ^ together with |lO| and |ll| we can immediately 
extract the scaling parameter A and the opening angle in 
terms of the contact parameter a and the shape param- 
eter m 



A(cr, m) = R 



sn (cr I to) dn (cr | m) 



a (cr, to) 
± 



2E{a\m)-a 
arccos [± (2dn^ (cr I to) - 1) 
— sign (2E (cr | to) — cr) 



(13) 

(14) 
(15) 



^ Some useful formulas and relations for the eliptic functions and 
integrals are briefly sketched in [ |25| and found in pq] in full 
depth. 

* Usually the parameter m is artificially assumed to be confined 
to < m < 1 but by the Jacobi's real transform for elliptic 
functions |26| they stay well-defined even for m > 1. 

^ A more visual parameter set [a, m, A) using the opening angle 
a = a{a, m) produces technical problems with non-uniqueness 
of loop representation. 



Plugging this into Eq. || we obtain the final form of the 
implicit constraint 



AL 

Jr 



sn (cr I to) dn (cr | to) 



2^; (cr I to) - cr 
arccos [± {2dn^ (cr I to) - 1)] 



(16) 



The curvature k (s) and the bending energy Eq. ^ follow 
from the explicit solution Eq. ^ to be 



k{s) 

^bend 



A 

AmA 

A 
4A 
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cn^ {t I m) dt 



— [(m- l)cr + £:(cr I m) 
A 



(17) 
(18) 
(19) 



The latter expression together with Eqs. [l], § - gives 
a lengthy expressio n fo r the total energy with the sign 



chosen ± as in Eq. 11 A 

Utoticr,m) = 



AA 



[2E (cr I to) - cr] [E (fJ I to) -I- (to, - 1) cr] 



sn (cr I to) dn (cr | to) 
+2Reads arccos [± (2dn^ (cr | to) - l)] (20) 

Now our problem of finding the ground state loop for 
given excess length AL reduces to a two variable (cr, to) 
minimization of Eq. |2^ under the constraint Eq. This 
final step is easily performed numerically. 



III. LOOP ZOOLOGY: SIMPLE AND CROSSED 
LOOPS 

We can scan now through the cr — to parameter plane 
and look at the shapes of the solutions and their en- 
ergies. In Fig. 3 we see a small (but most important) 
part of the whole parameter space and the correspond- 
ing different loop geometries. The dashed lines indicate 
parameter values which lead to constant excess length 
AL = 10 X 3.4nTO (corresponding to 100 bps) in accor- 
dance with the constraint Eq. |l6|. The shapes 1-7 are 
examples of lOObp-loops with different geometries. The 
whole parameter plane is subdivided by separation lines 
(solid) into regions of structurally different solutions. 
The large region starting at cr = contains exclusively 
simple loops (like 1,2 and 3) without self-intersections 
and nucleosome penetration. Above that simple-loop- 
region we find loops with a single self-intersection (4,5,6) 
and to the right the loops penetrate the nucleosome, like 
loop 10. There are also three other regions with single 
and double crossing points (7,8,9) where the loop can also 
be on the "wrong" side of the nucleosome like in 7 and 
8. 

We are interested in the energy minimizing loops and 
the underlying minimal energies as functions of the excess 
length AL. A density plot of these energies as function 
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FIG. 3: The set of possible ground-state solutions is charac- 
terized by two parameters, the contact point parameter cr and 
the loop shape parameter m. Solutions with constant excess 
length AL (here 10x3.4nm) are located along the dashed lines 
(e.g. loops 1-7). The solid lines separate loops with different 
geometric characteristics: simple (1,2,3), crossed (4,5,6) and 
"exotic" (7,8,9,10) loop shapes. 



FIG. 5: The ground state loop energy plotted vs the excess 
length AL. Note the energy maximum occuring for shorter 
loops. For much longer loops (around AL = 60nm) a tran- 
sition from simple uncrossed to crossed loop shapes occurs 
leading to a kink in C/min (AL). In the regime of low AL < Ip 
the elastic energy prevails strongly over entropy whereas for 
large loops the entropy starts to dominate the behavior pro- 
ducing a shallow energy minimum in the cross-over regime 
which roughly defines the predominant loop size. 
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o (contact param.) 

FIG. 4: Density plot of the total loop energy Eq. ^ (grayscale 
levelsets) as a function of a and m (same parameter range 
as in Fig. 3). The white contours denote lines of con- 
stant excess length AL = 1,2,3,5,10,20,50 x 3.4 nm. For 
given excess length the ground state is the point on the 
corresponding white line with the darkest background (note 
the different branches for given AL). The parameters are 
£ads ~ Q-lksT /nm and A = 50nm x ksT and R — 4nm. 



of the parameters a and m together with the correspond- 
ing lines of constant AL (with AL — 1, 2, .., 50 x 3.4nm) 
is given in Fig. 4. As can be seen from Fig. 3 there 
are, for a given AL, different branches of (cr, m) values 
corresponding to uncrossed, simply crossed and other ex- 
otic structures. Of all these structures for short excess 
lengths, AL ^ 20 x 3.4nm, the energetically dominant 



ones are simple (uncrossed) loops which we study first. 
Loops with larger excess length form crossed structures 
and are studied in Section 3.2. 



A. Simple Loops 

For simple uncrossed loops it is a straightforward nu- 
merical task to minimize Eq. 20 under the constraint of 
constant excess length, Eq. |16|. For Sads = 0-7kBT/nm 
and all the other parameters as above {A — 5Qnm x fc^T, 
R — 4nm) the ground state energy Umin as a func- 
tion^ of the excess length AL is shown in Fig. 5 (for 
AL ^ 60nrn; for longer AL- values crossed loops are more 
favorable as discussed in the next section). Remarkably 
we find that the loop energy is non- monotonous: For 

1/3 

small AL Umin increases with AL as (AL) ' (in accor- 
dance with Ref. [p^ where only small loops were stud- 
ied). At some critical excess length AL = ALcrit (which 
is approximately ALcrit ~ 2.2 x 3.4nm for Sads = 0.7 
ksT/nm) the loop energy reaches a maximum (here 
Umin{ ALcrit) ~ 26fcsT). Bcyoud that the energy de- 
creases with increasing AL. 

In the following we show how this behavior can be 
explained on the basis of the loop geometry. Naively 



^ Formally the quantisation condition AL = 1, 2, ... X SAnm holds 
as mentioned above. Nevertheless for clarity we consider the 
values in between as well. 
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FIG. 6; Two generic types of simple loop geometries (in the 
circle-line approximation): a) the subcritical loop with open- 
ing angle a < 7r/2 and b) the supercritical loop with a > 7r/2. 
In the former case the introduction of further excess length 
leads to an energy increase but in the latter case to a re- 
laxation of stress: The introduction of additional length at 
points Xl and Xr followed by a relaxation of the structure 
obviously decreases the total energy. 

one might argue as follows: For excess lengths shorter 
than the persistence length of DNA it is increasingly dif- 
ficult to store additional length into the loop because 
it requires increasing DNA deformation. On the other 
hand, for loops longer than Ip the bending energy con- 
tribution becomes very small and hence one expects such 
ground state loops relaxing with increasing AL. However 
the reason for occurence of a maximum of Umin around 
K, 2 excess DNA lengths, a value which is considerably 
smaller than the persistence length, is not obvious. In or- 
der to understand this finding one has to go beyond the 
simple handwaving heuristics and needs to take a close 
look at the details of the loop geometry. 

To this end we introduce here a simple approximation 
technique which leads to explicit expressions which can 
be more easily handled than the exact yet complicated 
expressions given above. We call this method the circle- 
line approximation and give a detailed exposition in the 
Appendix. As we will see this method is quite accurate 
and at the same time very intuitive. 

Looking at the geometrical shapes of the loops in Fig. 3 
we notice that each of them is subdivided into several 
sections of very high and very low curvature (cf. also 
Eq.|l^). In first approximation we replace the high cur- 
vature regions by sections of circles, the low curvature 
regions by straight lines (cf. Fig. 6). Furthermore, to 
keep the smoothness we assume that the lines are tan- 
gents to the circles. Generally in order to have reasonable 
approximations of all possible loop shapes we would need 
to consider compositions of several circles and lines (cf. 
for instance loops 3, 6, 7). However, if the adsorption en- 
ergies are not to high, i.e., if the opening angle a is "soft 
enough" and docs not impose such a severe bending like 
in loop 3, such multiply bent loops will not be relevant as 
ground state solutions. As it turns out for our problem 
we already obtain a quite good approximation by assum- 
ing that the loop consists of a single circular arc and two 
lines only. It is characterized by two quantities : 1) the 
arc radius r and 2) the nucleosome opening angle a (cf. 



Fig. 6 and Appendix). With these assumptions and after 
some elementary geometry the constraint Eq.|^ becomes 
simply 

AL = 2 (i? - r) (tan a-a) = const. (21) 

Note that the (more complex) second constraint Eq. || is 
eliminated through the "ansatz" per se. The total loop 
energy is given in terms of the loop radius r and the 
opening angle a 

ct 

Utot {a, r) = A V 2aR£ads 

r 

and by applying the constraint Eq.|l] (which this time 
can be solved explicitly!) we obtain Utot in terms of a and 
given AL 

, ^ „ tana — a „ \ 

- 2" (-^ 2i^(tana-a)-AL + ^^^^ ) ^''^ 

which is explicit in a. We note that this approximation 
for Utot is only reasonable for 2R (tana — a) > AL, i.e., 
for not too small a (vs. AL), otherwise the bending con- 
tribution diverges or becomes even negative (the latter 
is obviously absurd). The reason for this is that for very 
small angles a (compared to AL) uncrossed^ circle-line 
loops cannot exist for geometrical reasons. There this 
most basic approximation breaks down and we would 
have to approximate the loop by more than one circular 
segment. But as mentioned above, such loops (a small 
compared to AL) are not candidates for the ground state 
for moderate Eads ~ 0(1), and we therefore dispense 
with giving a discussion of this case. 

The nice thing about Eq. |2^ is that despite its simplic- 
ity and approximate nature it reproduces the position of 
the maximum in Fig. 5 quite well. We find the condi- 
tion for the critical excess length ALcrit from a simple 
geometric distinction between two loop shapes: the sub- 
critical loop (Fig. 6a) with its tangents not being parallel 
to the Y axis (a = 0) and the supercritical loop (Fig. 6b) 
having two or more tangents parallel to the line a = 0. 
Suppose now we add excess length to a subcritical loop by 
keeping the angle a — const. Obviously the loop-energy 
increases because the loop radius r becomes smaller. On 
the other hand in the supercritical case we have the oppo- 
site situation: the loop energy decreases with increasing 
AL. This is simply because we could cut the loop at two 
points {Xl and Xp in Fig. 6), introduce there the ad- 
ditional length (without changing the energy) and then 
relax the shape by letting it evolve to the new equilibrium 
while keeping a — const. Thus we can obtain the con- 
dition for the critical excess length ALcrit by assuming 
that the corresponding minimum a""" of Utot just crosses 



^ In contrary for crossed loops there still are solutions for small a 
(cf. next chapter). 
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the critical line 7r/2 line, i.e., a" 
searched ALcrit- 



(AL^"*) ^ 7r/2forthe 



d 

da 



Utot (a) = 



(23) 



Q=ir/2 



which can be solved for AL"^ 



4i? 8i?3 



(24) 



The latter can now be inserted in Eq. leading to 



Tjcrit 



■kA 
2R 



T^RSads 



(25) 



For the given values of R, A, Sads {R = ^nm , A = 50 
nrnkeT, Sads = 0.7 kBT/nm) we obtain AL'^"* = 7. 
37nni and Ui^l^ — 28. AksT which is in satisfactory 
agreement with the exact numeric results (AL'^"* = 7. 
19nm, Ulr^l'^ = 26. TksT). More generally, for not to 
high adsorption energies {£ads — 0.5 — 2.0 kBT/nm) the 
circle- line approximation works well and Eqs. 24 and |2^ 
reproduce the exact positions of the critical point typi- 
cally with a 5-15% accuracy. 

For an explicit parametric representation of the min- 
imal energy curve within the circle-line approximation, 
which in particular implies the upper results, the reader 
is referred to the appendix where the usefulness of this 
approach is also demonstrated for some other examples. 



B. Crossed and Entropic Loops 

A closer inspection of Fig. 4 shows that the ground 
state of loops switches from simple uncrossed loops to 
crossed loops when one reaches an excess length around 
50 nm. However, as can be seen for the crossed structures 
4, 5 and 6 in Fig. 3 these loops have a self-penetration 
at the crossing point. Therefore, a planar theory is in 
principle not sufficient to describe such structures. One 
possible formal cure for this problem would be to leave 
the plane and to consider the rod's self-contacts with the 
corresponding point-forces etc. in 3D as done by Cole- 
man et al. in a general theory of rod self-contacts [ p7| . 
However such a procedure leads to a significant loss of 
transparence, not only because of the third dimension 
entering the scene but also due to the necessity to sub- 
divide the rod into different regions with different forces 
acting in each of them. Instead of following Coleman at 
al. p7| we decided to treat the self-interaction in a per- 
turbative manner as follows. If the self-contact point is 
not too close to the nucleosome the rod is not severely 
deflected out of the plane by its self-interaction. Thus it 
remains roughly planar with some out of plane bending in 
Z-direction of the rod sections between the nucleosome 
and the crossing point. This will cost some additional 
bending energy Udef that is roughly given by (cf. Ap- 
pendix) 



Udef (c, m) 



, IP tan a.((T,m) 

R tan^ a (u, m) — 



(26) 



Here p := d/R with d k, Inm is the DNA radius. We 
neglect the slight twisting of the rod induced by the non- 
planarity of the DNA and consider the bending only. The 
deflection energy Eq. |2^ can be phenomenologically in- 
corporated into the model by simply adding it to Eq. ^ 
as a correction term to obtain the final form of the total 
energy U^^^ 



Utot (c^' ™) 



Utot (f, to) for uncrossed (simple) loops 
Utot (c, m) + Udef (c, m) for crossed loops 



With this additional modification of Utot we computed 
numerically the minimal energy (ground state) solution 
for any given excess length AL. The graph of the ground 
state energy versus AL is shown if Fig. 5. We find 
that even with the inclusion of the out-of-plane deflec- 
tion there is still a critical length ALcross (here « 60nm) 
where the crossed loops become energetically more fa- 
vorable than the simple uncrossed. This behavior that 
we call the "crossing transition" can be rationalized by 
noting that for long enough loops the adsorption energy 
(proportional to a) starts to dominate over the bending 
energy so that loops with smaller a become increasingly 
favorable. From the critical length ALcross on, the gain 
in adsorption energy (by diminishing a) is more than suf- 
ficient to outweigh the (slight) increase in bending energy 
together with the additional self-interaction term, Eq.|2q. 

Increasing the length even further we leave the elas- 
tic energy dominated regime in which the entropic ef- 
fects can be neglected due to short loop length (< persis- 
tence length). For larger lengths entropic effects become 
more and more important and we ultimatively enter the 
entropic loop regime. The crossover between these two 
regimes is hard to handle analytically ; for the case of 
closed loops a perturbative description has been given in 
Ref. . For our purpose it is sufhcient only to consider 
the asymptotic behavior. In the large loop limit where 
the loop is longer than several I p the chain looses its " ori- 
entational memory" exponentially and behaves roughly 
as a random walk which starts from and returns to the 
same point. The entropic cost for gluing the ends of a 
random walk (long loop) together is then given by 



U = 3/2kBT\n{AL/lp) + Eq + So 



(27) 



The first constant, Eq w 6.5 kpT is the bending -I- ad- 
sorption energy contribution of the overcrossing DNA 
segments leaving / entering the nucleosome which can be 
determined by minimizing the crossed loop energy (cf. 
Appendix Eq. ^) for AL oo. The second additive 
constant 5*0 ~ 0{kBT) accounts for the entropic contri- 
bution of DNA-histone octamer interaction volume (the 
proximity necessary for the histone octamer and DNA to 
see each other) . Although the latter constant is not easy 
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to estimate the following prediction is not sensitive to any 
additive constant. We expect a free energy minimum to 
occur at the overlap between the elastic (AL < Ip) and 
entropic (AL ^ Ip) region where the decreasing elastic 
energy is overtaken by the increasing entropic contribu- 
tion. 

The free energy, Eq. leads to an algebraically de- 
caying probability w (Ai) for the jump lengths scaling 

as w ^ (Ai)^^^^. In general, power law distributions 
of the form w ~ (AL)^^ with 7 > 1 lead to superdif- 
fusive behavior of the random walker (here the nucle- 
osome). According to Levy's limit theorem the proba- 
bility distribution of the random walker (more precisely, 
the distribution of the sums of independent random vari- 
able drawn out from the same probability distribution 
w ~ (AL)"'') converges to a stable Levy distribution of 
index 7 — 1 |3^, |3|, |3|] . This so-called Levy- flight dif- 
fers in many respects from the usual diffusion process as 
for short time intervals big jumps are still available with 
significant probability. Moreover, all moments (besides 
possibly the first few ones) diverge. For our case 7 = 3/2 
even the first moment does not exist. We note that the 
value 3/2 is based on the assumption of an ideal chain (no 
excluded volume); in general the excluded volume leads 
to self-avoiding-walk statistics with a slighly larger value 
of 7 around 2.2 |^ (cf. also Ref. |3^). In that case one 
has a finite value of the first moment, i.e., of the average 
jump length. 

IV. THE DYNAMICS OF NUCLEOSOME 
REPOSITIONING 

In the preceding sections we have computed the typi- 
cal energies involved in the formation of arbitrary sized 
loops. Assuming that a slow creation followed by a fast 
termal migration of loops around the nucleosome is the 
governing mechanism for nucleosome repositioning we 
start now considering the repositioning dynamics. In or- 
der to describe the time-dependent evolution of the nu- 
cleosome position we consider its probability distribution 
along a DNA segment of a length N x IQbp and write the 
master equation governing the jump process 

d ^ iv 

—Pi = Wj.pj - Pi Wij (28) 

where pi is the probability for the nucleosome being at 
the admissible® position i on the DNA segment. The 
transition rate matrix W_ = {wij ) is given by 

_ j Ca exp (-j^U^in (ho \i - for i ^ j 

(29) 



Spaced by a multiple of 10 bp from the initial position 




FIG. 7: Relaxation dynamics of two initial states of nucleo- 
some positions on a short DNA segment (147 -f 90 bp): a) 
the nucleosome starting from an end and b) the nucleosome 
starting from the middle position. The time unit is the inverse 
Arrhenius activation factor C~^^ (compare text). 

where hpt = iAnm (DNA helical pitch). Ca denotes the 
Arrhenius constant involved in the loop formation pro- 
cess that has in principle to be determined experimen- 
tally. The rough estimate of C^^ — lO^^s is provided 
in Ref. where it was shown that Ca is essentially 
given by the inverse lifetime of the loop (denoted by A in 
that paper) . This means that typical repositioning times 
range from seconds to hours. 

The (formal) explicit solution of Eqs. ^ together 
with the previously obtained minimal energy C/min is 
given by 

p (t) = exp(^t)p (0) 

The latter solution can now be considered for different 
cases: for short or long DNA chains and for the nucle- 
osome placed in the middle or at the end of the chain. 

For short DNA segments we expect a slow reposition- 
ing rate due to high energies involved in small loop for- 
mation. In Fig. 7 we depict the repositioning of a nu- 
cleosome on a DNA piece of a length 147-1-90 bp. Start- 
ing from an end positioned nucleosome (Fig. 7a) we ob- 
serve a behavior that is completely unlike a local dif- 
fusion mechanism: the jumps bigger than w 2 x SAnm 
start to dominate over the smaller local ones, which fol- 
lows from the loop formation energy cf. Fig. 5. Con- 
sequently, in the initial phase of repositioning (of such 
an end-positioned population) the nucleosomes will pre- 
dominantly jump between the two end positions. Later, 
on a much larger timescale they gradually start to ex- 
plore the positions towards the middle of the DNA seg- 
ment. Could we extract such a behavior from an experi- 
ment using gel-electrophoretic separation (as in j^, [^)? 
The basis of such separations is the fact that the gel- 
electrophoretic mobility of nucleosomes on DNA pieces 
(longer than 147bp) increases roughly linearly with its 
distance from the middle position, i.e., DNA pieces with 
the nucleosome sitting close to the end run much faster in 
gels than equivalent middle positioned nucleosomes do. 
We can exploit this (empirical) fact to mimic the outcome 
of a gel-electrophoresis experiment (cf. Figs. 8 and 10). 
In Fig. 8a we depict such a simulated gel pattern for the 
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1 2 3 4 5 




1 2 3 4 5 

FIG. 8: Typical (1-D) gel electrophoresis signatures expected 
for the relaxation dynamics of the two species from Fig. 7: 
a) nucleosome starts from an end and b) from the mid- 
dle position. The lanes 1-5 correspond to incubation times 
(l,5,10,20,100)xl0* CJ^ respectively. Note: the population 
of distant bands in b) lanes 2-4 occurs first, in sharp contrast 
to what we expect from a simple (local) diffusive behavior. 



middle positioned nucleosome. Since symmetric species 
are not distinguished by this experimental method and 
are projected onto the same bands (symmetric left/right 
positions lead to the same mobility) , the expected nonlo- 
cality of motion cannot be extracted from the structure 
of the bands. 

For the same short segment, but with the nucleosome 
starting from the middle position (Fig. 7b) the situation 
is slightly different: the neighboring positions are pop- 
ulated more homogeneously, although there is a small 
initial underpopulation of the 2 x 3.4nm distant position 
as expected from the energy maximum occurring there. 
In this case, a slight initial "population gap" can be ob- 
served in gel electrophoresis (Fig. 8b) which in this case 
would be sufficient to distinguish between a jumpy and a 
diffusive behavior, since the latter would obviously lack 
the " population gap" . 

In the case of longer DNA (but still not entropic seg- 
ments) like the 147-1-300 bp segment in Figs. 9 and 10, 
similar effects as for the short segments are expected but 
with significantly faster relaxation times by typically 2-3 
orders of magnitude as compared to the corresponding 
short segment populations. The corresponding (simu- 
lated) electrophoretic gels are shown in Fig. 10 where for 
the centrally positioned case (Fig 10b) the "population 
gap" effect is even more pronounced than in the short 
segment case. 

For even longer DNA segments we expect the gap effect 
to persist (data not shown) and the optimal jump size to 
be around 2-3 x I p corresponding to the free energy min- 
imum in Fig. 5. For very long DNA segments, the nu- 
cleosome repositioning behavior implied by the big-loop- 



mechanism becomes strongly non-local which contrasts 
a local diffusive motion as expected from cork-screwing 
motion (cf. Refs. Q |^, ^1 ) or small loop reposition- 
ing as considered by Ref . llOl . As mentioned above, this 
superdiffusive behavior has diverging moments which im- 
plies strongly enhanced nucleosome transport along very 
long DNA pieces. However such an ideal superdiffusion 
of nucleosomes could hardly occur in vivo because free 
DNA segments between subsequent nucleosomes (DNA 
linkers) are never longer than ~ 0(lp). Furthermore 
the neighboring nucleosomes might be a significant bar- 
riers (if not for loop formation then) for loop migration 
around the nucleosome which is an indispensable event 
for loop-mediated repositioning. 



V. CONCLUSIONS AND DISCUSSION 

In this study we examined a possible mechanism for the 
repositioning of nucleosomes along DNA which is based 
on the formation and diffusion of intranuclcosomal loops. 
The most important outcome of this study is the predic- 
tion of two classes of loops that might occur: (1) small 
lObp- loops and (2) large loops with a wide distribution 
of stored lengths with a weak peak at roughly two times 
the DNA persistence length. 

The small loops were already discussed in Ref. [0 and 
led to the prediction of repositioning steps of 10bps. Fur- 
thermore, the repositioning time should be of the order 
of an hour, a consequence of the large activation energy 
required to form a loop. This might explain the strong 
temperature dependence of the typical repositioning time 
1^. In fact, by lowering the temperature from 37° to 4°C 
no redistribution within one hour was detected in that ex- 
periments. Assuming a loop formation energy of 23kBT 
one finds indeed a slowing down of this process by factor 
of 13. 

On the other hand, the large loop repositioning consid- 
ered here turns out to be energetically much more favor- 
able. Loops with an extra length of 21 p have an energy 
that is roughly 12-13A:bT smaller than that of a lObp- 
loop. To a certain extend this is because such loops can 
have a very small nucleosome opening angle by forming 
crossed loops but the main contribution stems from the 
significantly decreased DNA bending energy. One there- 
fore expects that repositioning via large loops should be 
the dominant process on sufficiently large DNA pieces 
and that the typical times are much shorter than the one 
for small loop repositioning (say, of the order of minutes) . 

So far, however, the experiments did not report such 
events. Meerseman et al. ^, for instance, found on 
short DNA pieces of 207bps length results that are con- 
sistent with lObp repositioning ~ as we would expect for 
such short DNA fragments. However, when they redid 
the experiment with a 414bp long piece, a tandem repeat 
of the 207bp DNA, their analysis of the complicated band 
patterns observed in 2D gel electrophoresis did not show 
any indication that the nucleosome was able to move from 
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FIG. 9: Relaxation dynamics of two initial states of nucleosome positions on a longer DNA segment (147 + 300 bp): a) end 
positioned and b) centrally positioned initial species. Note the initial difference in relaxation timescales for a) and b) (which 
are due to different loop energies involved). 




FIG. 10: The (1-D) gel electrophoresis signatures simulated for the relaxation dynamics of the two initial species from Fig. 9. a) 
End positioned (lanes 1-5 corresponding to incubation times (1,2,3,10,50) x lO** C^^ ) and b) centrally positioned (incubation 
times (1,2,3,10,50) X 10** CJ^). 



one half to the other. 

Hence, the question arises if the repositioning observed 
in these experiments was facilitated via the loop mecha- 
nism or if it occurred via a different process. An analysis 
of the results is made especially difficult by two com- 
plications: (a) the nucleosomes seem to prefer to sit on 
the ends of the DNA fragments and (b) most of the ex- 
periments use strong positioning sequences (like the 5S 
rDNA sequence). This means that, independent of what 
the repositioning process might be, the nucleosomes have 
certain preferred positions and these might obscure the 
underlying repositioning process. 



With regard to this fact, let us consider two other repo- 
sitioning mechanisms that one could imagine. The first 
one is that the nucleosome detaches completely from the 
DNA and attaches at some other position (or even a dif- 
ferent DNA molecule). This process, however, seems to 
be excluded by two facts (among others). First that no 
repositioning from one half to the other of the 414bp 
DNA or to competitor DNA fragments was observed 

1^ . Secondly, once completely detached from the DNA 
template the histone octamer becomes unstable and dis- 
integrates into a tetrameric and two dimeric subunits 
which makes an effective nucleosome reconstitution difh- 
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cult. 

The other mechanism could be a local screwing motion 
as already suggested in Ref . . This process would lead 
to a repositioning with one bp per step. The preponder- 
ance of lObp steps observed for the 5S rDNA experiments 
could then be explained as being due to the fact that the 
positioning sequence prefers the nucleosome rotationally 
positioned on one side of the DNA where it can be easily 
bent around the octamer. Also 10bps (and even a few 
multiples of 10bps) apart this effect can still be seen and 
hence the nucleosome would prefer positions multiples of 
10bps apart. To our best knowledge, the experiments to 
date do not allow to distinguish whether the lObp repo- 
sitioning works via small loops or via cork-screwing. 

It would be therefore important to perform experi- 
ments on DNA pieces that do not provide the nucleo- 
some with a preferred rotational setting. In that case 
the lObp footprint should disappear if nucleosomes repo- 
sition themselves via cork-screwing. It would also be im- 
portant to perform experiments with rather long DNA 
fragments since we expect that large-loop repositioning 
can be detected in such systems. 

Finally, we note that nucleosome repositioning in vivo 
is facilitated via so-called chromatin remodeling com- 
plexes, huge multi-protein complexes that harness energy 
by burning ATP ||, |, |. There are basically two ma- 
jor classes: ISWI and SWI/SNF. The first one seems to 
induce small scale repositioning which might work via 
twisting DNA that leads to a corkscrew movement as 
discussed above. It might, however, also be possible that 
this complex induces small loops on the nucleosome as 
recent experiments on nicked DNA suggest [Q. The 
other class of remodeling complexes seems to induce large 
loop structures as they have been observed recently via 
electron spectroscopy |35| . Whatever the details of the 
functions of these remodeling complexes might be, it is 
tempting to speculate that they catalyze and direct pro- 
cesses which might even take place when they are not 
present - like small loop and large loop formation as well 
as screwing. In this case the computed looping energy 
(cf. Fig. 5) and repositioning rates might give a first hint 
about ATP requirements and the dynamics of enzymatic 
repositioning. 

Another interesting and very prominent system known 
to mediate nucleosome repositioning via loop formation 
is unexpectedly the ubiquitous RNA-Polymerase (RNA- 
P). It is found to be able to transcribe DNA through nu- 
cleosomes without disrupting their structure, yet moving 
them upstream the DNA template, i.e., in the opposite 
direction of transcription Jsof . To rationalize this seem- 
ingly paradoxical finding Felsenfeld et al. introduced a 
DNA looping model which assumes that the RNA-P 
crosses the nucleosome in a loop. This would indeed ex- 
plain the backwards directionality of repositioning. An 
interesting question in this context is how our intranucle- 
osomal loops considered above relate to those formed by 
the RNA-P. Can we say something about the reposition- 
ing distance distribution, does the looping energy (Fig. 5) 



apply here? The geometry of RNA-P - DNA complex on 
a nucleosome is certainly different from the simple loop 
case, as ingoing and outgoing DNA from RNA-P enclose 
a (rather soft yet) preferential angle of « 100° (depen- 
dent on RNA-P type, cf. Refs. ||, H ||). The latter 
facilitates the loop formation as the free DNA has to 
bend less to fold back onto the octamer surface. Besides 
the apparent differences from the " naked" intranucleoso- 
mal loops problem, a slight generalization of our present 
model which incorporates the preferential RNA-P "open- 
ing" angle can be performed within the same mathemat- 
ical framework developed here. It would be interesting 
to compute the resulting nucleosome transfer distance 
on short and long DNA templates in an analogous man- 
ner as performed above. An outcome of such a study 
could be, for instance, an answer to questions like: what 
is the highest linear nucleosomal density in polynucle- 
osomal arrays, up to which nucleosomes are not to be 
removed from the DNA template (due to loop formation 
and nucleosome transfer prohibited by the neighboring 
nucleosome) during transcription. 

Such fundamental biological questions make a further 
elaboration of intra-nucleosomal loop theory, its gener- 
alization to different loop geometries, and finally its ap- 
plication to different loop creating proteins (SWI/SNF, 
RNA-P) an intriguing task for future work. 



VI. APPENDIX: THE CIRCLE-LINE 
APPROXIMATION 

Although the Kirchhoff 's analogy provides us with es- 
sentially analytic solutions for the rod deformed in plane, 
the occurrence of boundary conditions (like Eqs. || and 
prevents us in most cases from obtaining analytical 
expressions of all the parameters characterizing the solu- 
tion (like cr and m above) . To overcome this problem, we 
suggest here a simple geometric approximation scheme 
which will prove to be useful in obtaining analytic re- 
sults for loops within a reasonable accuracy (usually with 
a deviation of 5-15% from the exact numeric results). 

The main idea is the following. The curvature and 
the energy (Eqs. ^ and |l^) of the loop contains the 
cn(cr|r7i) function which for < to < 1 has the typi- 
cal oscillatory behavior depicted in Fig. 11 (left). This 
suggests to approximate the curvature function simply 
by a step function consisting of an alternating sequence 
of negative, zero and positive piecewise constant curva- 
tures. Consequently the corresponding rod shape (Fig. 11 
right) is approximated by a sequence of circles (positive / 
negative constant curvature) and lines (zero curvature). 
An analogous approximation procedure can also be per- 
formed in the case to > 1 where the cn function has a 
natural analytical continuation into a dn function with a 
modified second argument (cf. Ref. [p6[). 

Using this approximation ansatz several problems con- 
cerning planar rods reduce to elementary geometry as 
seen from the following simple but illustrative examples. 
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FIG. 11: The circle-line approximation for planar rods. 
The curvature of an equihbrium rod shape (cn-function, cf. 
Eq. ^) is approximated by a periodic sequence of step- 
functions. The latter corresponds to an approximation of the 
rod shape by a sequence of straight lines (k = 0) and circles 
(k =const.) glued together in a smooth manner (continuous 
tangents). 





FIG. 12: Three applications of the circle-line approxima- 
tion: Problems with complex constraints reduce to simple ge- 
ometries leading to good approximations: a) the Yamakawa- 
Stockmayer angle b) simple loops and c) crossed loops (see 
the Appendix text for details) . 




1) The Yamakawa-Stockmayer angle |29 : Two points 
on the rod are glued together without restricting the ori- 
entation of the tangents, e.g., a protein connects two dis- 
tant points on DNA (cf. Fig. 12a). What is the pre- 
ferred angle x between the tangents in the ground state 
of the rod? By imposing a fixed total rod length L we 
have the simple constraint L — (2 cot f + X + tt) r from 
which we can eliminate r and write the elastic energy of 
the configuration as U^^\ — ^ (x + ^) (2 cot ^ -I- x -I- tt) . 
Its minimization leads to the transcendent condition 
Xmin + TT = tan Xmin wlth the only relevant solution 
Xmin ~ 77.5°. The latter angle differs by 5% from the 
exact result Xmin ~ 81.6° (by Yamakawa and Stockmayer 
in 1^) which is satisfactory regarding the simplicity of 
the computation. 

2) Simple and crossed loops (Fig. 12 b,c): We can eas- 
ily derive an approximate energy expression for simple / 
crossed loops as a function of the excess length Ai and 
the opening angle a. By applying simple geometry the 
excess length constraint can be easily eliminated (the tan- 
gency constraint is trivially fulfilled by the ansatz) and 
we arrive at 



Ustrap [a) = 2a[A 



tan a — a 



2R{i&\ia -a)- AL 



RSa. 



ds 



for simple loops and 



(30) 



(a) 



2a A 



tan a — a 



RSad 



A L [nm] 



AL - 2i?(tana - a) 
+Udef (a) (31) 

for crossed loops where A, R and Eads defined as above 
and Udef being the excluded volume interaction at the 
crossing point, which is considered below (and applied 
in the main text as Eq. We remark that the above 
expressions for Usimp and Ucross are valid within certain 
a intervals which are given by the restriction < a < tt 
and by the condition that the first terms in the brackets 
of Eqs. ^ and ^ are positive (these are the necessarily 
positive bending energy contributions in the two cases.) 

These fairly simple expressions can now be used in the 
two cases to obtain explicitly the ground state energies 
by minimizing Eq. |3^ and Eq. ^ with respect to a. For 
instance, setting U^^^^ [a) = we obtain a transcenden- 
tal equation for a. We can now use the fact that this 
condition is algebraic in AL so that we can solve it for 
AL = AL (a). Thus instead of finding a = a{L) (which 
cannot be given in an explicit form) we obtain explicitely 
its inverse: 

AL (a) _{2~c)G [a) + cH {a) 



R 



1 



FIG. 13: Comparison of the adsorption and bending energy 
contributions {Uada and Ubend) as well as the total ground 
state energy Utot of the simple loop. The fat lines represent 
the circle-line approximation (cf. Eq. ^) whereas the thin 
lines show the corresponding exact expressions, Eqs. ^ and 
po| (thin line). The parameters are Eads ~ O.TfcsT/nm and 
A = 50nm x ksT and R = 4nm. 



X 



^ [(2 -c)G {a) + cH {a)f - 4(1 - c)G2 (a) 



1 -c 

with the abbreviations 

G {a) = tana — a 
H (a) = atan^ (a) 



(32) 
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FIG. 14: The out of plane deflection of the incoming/outgoing 
DNA due to excluded volume in the top projection (left) and 
seen from the side (right). In the latter case (for the sake of 
visual clarity) the two rods are depicted in a single plane, i.e., 
rotated around their contact point (grey dot). 



In Eq. |3^ the introduced dimensionless constant is c = 
(1 + 2R^eads/Ay^ (0 < c < 1, and c = 0.69 here) and 
X is the sign accounting for different branches of the a 
parametrized solution 



X 



-1 for < a < 7r/2 

±1 for Tt/2 <a< Omax (c) 



(33) 



Note that for a < 7r/2 there is only one branch but for 
a > 7t/2 we have two branches^ (il) for AL{a). The 
maximal opening angle amax (c) is obtained by setting 
the discriminant (expression below the square root) in 
Eq. ^ equal to 0. 

From Eq. ^ together with Eq. ^ we obtain an ex- 
plicit parametric representation of the minimal energy 
curve for simple loops. A comparison of the approximate 
minimal energies (Eq. |3^ and Eq. ^o|) with the exact min- 
imal energy (cf. also Fig. 5 for AL < 60nm) is shown in 
Fig. 13. We find that the quantitative agreement is quite 
satisfactory taking the simplicity of our ansatz into ac- 
count. We note here that analogous computations as we 
have shown for simple loops can be performed for crossed 
loops as well. 

For AL — !■ wc find after an appropriate expan- 
sion of Usimp around a = that the ground state en- 

1/3 

ergy scales as Usimp ~ {^L/R) in agreement with 
Ref. 0. Further we obtain the excess length at which 
the loop ground state energy is maximal by setting 
dUsimp (a) /da\a=^/2 = 0. From this follows the critical 



within 2-15% with the exact numerical result for a wide 
range of adsorption energies with deviations becoming 
larger for adsorption energies above Sads — 2.0 kBT/nm 
(data not shown). 

3) The overcrossing potential for crossed loops (Fig. 
14): The outgoing DNA path is perturbed out of the 
plane due to the interaction with the ingoing DNA (and 
vice versa in a symmetrical manner) . Because of that our 
simple planar and phantom model (no self interaction) 
needs modifications. Instead of solving this (nonpla- 
nar) problem within the general theory of self-interacting 
deformed rods as in Ref. (which is a feasible but 
rather technical numerical task) we can treat the out-of- 
plane deformation perturbationally. The first assump- 
tion we make here is that the overall shape of the crossed 
loop does not deviate much from a planar configuration 
though the orientation of its (effective) plane might be 
slightly deflected from the nucleosomal plane. Conse- 
quently the small perturbation out of the plane and the 
deformation in plane essentially decuple into a sum of 
two energy contributions as in Eq.^. Again by simple 
geometry (cf. Fig 13), the second (out of the plane) term 
in can in first approximation be written as 



d arctan 



Udef (a) 



= <^ 2A- 



( 2dx{a) \ 



(a) — 
oo otherwise 



for X (a) > d 



length ALcrit as discussed in the main text (cf . Eq. 24 ) 
This simple approximate expression for ALcrit agrees 



(34) 

where d w Inm is the thickness of DNA and x (a) :~ 
Ri&na the length of the crossed segment. In our sim- 
ple approximation the self-interaction energy diverges for 
X — > d -I- as ^A{x — d)~^ (extreme deformation) and 
approaches zero for x — > oo as 4:Ad^x~^ (weak defor- 
mation) . 

We finally note that besides the above given examples 
it is possible to apply the circle-line approximation to 
several other standard problems of rod theory like the 
first and especially the higher order Euler buckling in- 
stabilities to obtain qualitatively the known results from 
buckling theory with very little effort. Thus the circle- 
line approximation when applied appropriately turns out 
to be very useful and generally allows computationally 
inexpensive qualitative and quantitative insights into the 
behavior of (planary) deformed rods. 
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